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Context. The lower limit to the distribution of orbital periods P for the current population of close-in exoplanets shows a distinctive 
discontinuity located at approximately one Jovian mass. Most smaller planets have orbital periods longer than P ~ 2.5 days, while 
higher masses are found down to P ~ 1 day. 

Aims. We analyze whether this observed mass-period distribution could be explained in terms of the combined effects of stellar tides 
and the interactions of planets with an inner cavity in the gaseous disk. 

Methods. We performed a series of hydrodynamical simulations of the evolution of single-planet systems in a gaseous disk with an 
inner cavity mimicking the inner boundary of the disk. The subsequent tidal evolution is analyzed assuming that orbital eccentricities 
are small and stellar tides are dominant. 

Results. We find that most of the close-in exoplanet population is consistent with an inner edge of the protoplanetary disk being 
located at approximately P >2 days for solar-type stars, in addition to orbital decay having been caused by stellar tides with a specific 
tidal parameter on the order of ^ 10'. The data is broadly consistent with planets more massive than one Jupiter mass undergoing 
type II migration, crossing the gap, and finally halting at the interior 2/1 mean-motion resonance with the disk edge. Smaller planets 
do not open a gap in the disk and remain trapped in the cavity edge. CoRoT-7b appears detached from the remaining exoplanet 
population, apparently requiring additional evolutionary effects to explain its current mass and semimajor axis. 



1. Introduction 

Close-in planets (semimajor axis a < 0.1 AU) constitute a spe- 
cial subset of the exoplanetary population. Since it is unclear 
whether in-situ formation occurs, the current orbital and physical 
characteristics of these planets provide important constraints on 
their past evolution and formation process. Several mechanisms 
have been proposed to explain the pile-up of hot planets with a 
three day orbital period, including a truncation of the gaseous 
disk by the star (Lin et al. 1996, Kuchner and Lecar 2002), plan- 
etary scattering combined with Kozai resonance and tidal circu- 
larization (Nagasawa et al. 2008), planetary evaporation (Davis 
and Wheatly 2009), and tidal interactions with the parent star 
(Jackson et al. 2009). 

In particular, Kuchner and Lecar (2002) suggested that a 
giant planet in circular orbit could halt its migration when its 
orbital period was half that of the inner edge of the disk. In 
this configuration, all the planet's circular Lindblad resonances 
would lie in the inner cavity (IC) and no further interchange of 
angular momentum would take place. Masset et al. (2006) per- 
formed a series of hydrodynamical simulations to follow the evo- 
lution of low-mass planets in disks including an IC. They found 
that all bodies migrated until reaching a point slightly exterior to 
the cavity edge, where they were effectively trapped in a stable 
configuration in almost circular orbits. Although this result ap- 
pears different from that predicted by Kuchner and Lecar (2002), 
each is valid, as we shall see, for a different range of planetary 
masses. 

The first reference to a possible correlation between mass 
and orbital period for close-in planets was proposed by Mazeh 
et al. (2005) for only six transiting bodies. They found that both 
parameters seemed to follow a linear law, with more massive 
bodies being located at smaller semimajor axes. Southworth et 
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Fig. 1. Distribution of orbital periods and planetary masses for 
close-in exoplanets with orbital period P < 12 days. Data 
from http://exoplanets.eu Black circles denote planets with both 
Doppler and transit data, while gray circles mark bodies with- 
out detected transits. Empty squares correspond to planets with 
retrograde orbits with respect to the stellar spin. 



al. (2007) and Davis and Wheatley (2009) extended the analysis 
to a larger transiting population, finding a similar result although 
with a much broader dispersion. They proposed that smaller 
planets closer to the star might have been lost because of evap- 
oration, similar to that currently ongoing at least in HD209458b 
(Vidal-Madjar et al. 2003) and WASP-17 (Anderson et al. 2010). 
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Jackson et al. (2009) also analyzed the distribution of close- 
in planets, this time focusing on the correlation between the 
semimajor axis of the planet and the age of its star. They 
found that the lower limit of the semimajor axis was lower for 
younger stars, which implies that tidal effects could be respon- 
sible. Exoplanets with very short orbital periods in older stars 
would have had enough time to be tidally disrupted, thus they 
would only be presently observable in relatively young systems. 

In this paper, we revisit the mass-period distribution, tak- 
ing advantage of the recent increase in the exoplanet population. 
Figure[T]shows the orbital periods P, as a function of the mass m, 
for the known population of exoplanets with P < \2 days (137 
planets). Black circles correspond to cases for which both tran- 
sits and Doppler data are available; bodies without detected tran- 
sits are shown in gray. Exoplanets in apparent retrograde mo- 
tion with respect to stellar rotation are identified by an empty 
square. These are WASP-8b (Queloz et al. 2010), WASP-17b 
(Anderson et al. 2010), WASP-33b (ColHer Cameron 2010), 
Hat-P-7b (Winn et al. 2009), WASP-2b (Triaud et al. 2010), and 
WASP-15b and WASP-17b (Triaud et al. 2010). Although it may 
be argued that these bodies are not consistent with planetary mi- 
gration (Triaud et al. 2010), they may also point towards primor- 
dial spin-orbit misalignment and not be related to subsequent 
orbital evolution of the planets (Lai et al. 2010). 

The distribution exhibits a noticeable "step", exoplanets 
larger than one Jupiter mass (Mjup) appear to have a lower inner 
boundary (down to ~ 1 day) while for m < Mjup the distribution 
seems restricted to larger values of P. The only exceptions are 
three bodies in the Super-Earth range, CoRoT-7b, GJ1214b, and 
GJ876d, which are all marked in the plot. Of these, the latter two 
planets belong to low-mass stars (M» - 0.17 and M» - 0.32 so- 
lar masses, respectively), thus constitute special cases. CoRoT- 
7b, however, belongs to a solar-type star (Rouan et al. 2009). 
This planet has a very short orbital period (~ 0.85 days) but also 
a very low mass {m ~ 0.015Mjup), and does not seem to comply 
with the rest of the exoplanet distribution. In particular, Jackson 
et al. (2009) pointed out that CoRoT-7b could reach the Roche 
radius on timescales of 10^ - lO** years, depending on the value 
of the specific tidal parameter Q'^. 

Regardless of these isolated cases, there seems to be a very 
clear discontinuity (or bump) in the mass-period distribution, lo- 
cated at approximately m - Mjup. Moreover, for high masses the 
lower limit in orbital periods appears very close to a 2/1 mean- 
motion resonance with the disk edge for small planetary bodies. 
This appears consistent with a scenario in which the planetary 
traps proposed by Masset et al. (2006) would dominate the low- 
mass region, while the mechanism of Kuchner and Lecar (2002) 
would be mainly responsible for the upper end of the mass spec- 
trum. 

The main objective of this study is to test whether the com- 
bined action of planetary traps in the gaseous disk plus subse- 
quent tidal interactions with the parent star could explain the 
observed distribution of close-in exoplanets. Since exoplanets in 
retrograde orbits should have exotic disk-planet and tidal evolu- 
tions, the study of these exosystems is beyond the scope of the 
present work, and we focus mainly on bodies believed to have 
orbital motion in the same direction as the stellar spin. Even in 
this case, we assume a zero incUnation with respect to the stellar 
equator. 

In Section 2, we present a series of hydrodynamical simu- 
lations adopting different planetary masses and analyzing the 
relative halting distance from the central star. Not only are we 
interested in seeing whether such a hybrid and mass-selective 
process is possible, but also whether the boundary between both 
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Fig. 2. Snapshot of the surface density profile of one of our sim- 
ulations. The inner cavity is centered around r = 1.8. In the 
top frame, the low-density regions are shown in black, while 
the high-density regions are shown in white. The planet (m = 
O.lMjiip) is located on the x-axis. 

mechanisms is consistent with the observed distribution of close- 
in planets. Section 3 is devoted to the subsequent evolution of 
exoplanets under the stellar tide and their effects on any initial 
disk-driven distribution in the mass-period diagram. In Section 
4, we analyze the case of the CoRoT-7 planetary system and 
present possible explanations of the present location of CoRoT- 
7b. Finally, conclusions close the paper in Section 5. 

2. Hydrodynamical simulations 

2. 1 . Initial conditions 

Our simulations were carried out using the FARGCQ hydro- 
code (Masset 2000a,b) with a central star of one solar mass. We 
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Fig. 3. Evolution of semimajor axis of two one-planet systems 
migrating and halting near the disk inner cavity edge (IC). Black 
curve corresponds to m - O.lMjup and gray to m - lOMjup. The 
IC is centered at r/c = 1.8, while its interior 2/1 mean-motion 
resonance (MMR) is located at r ^ 1.13. Time is in orbital peri- 
ods. 

consider two-dimensional, locally isothermal disks. The unit of 
length, ro, is arbitrary (our simulations are scale free and can 
be rescaled to disks of arbitrary sizes), but it can be considered, 
for our purpose, that ro ~ 0.01 AU. Our disks, prior to a ID 
relaxation dedicated to the creation of the IC, have a uniform 
surface density Z = 10"^ M^.r^^, which translates, in our case, 
to E ~ 900 g.cm"^. The asymptotic radius of our planets is virtu- 
ally independent of this value, which essentially prescribes how 
long it takes for them to reach it. In a similar manner, the kine- 
matic viscosity of our disks is uniform outside the cavity, and is 
set to V = IO^^TqQq', where Qq is the orbital frequency at ra- 
dius ro. Some initial test runs were performed with other values 
of 2 and v with no significant change in the results, except for 
the migration timescale. The disk aspect ratio is set uniformly to 
Hjr - 0.05, hence the value of the a parameter that character- 
izes turbulence (Shakura & Sunyaev, 1973) just outside the IC 
(see below) is a ~ 3 • lO"-'. 

Our disks are described on a polar grid with 384 zones in az- 
imuth and 306 zones in radius. The inner radius of the mesh is at 
0.2ro, while the outer radius is at 3.5ro. The IC was generated us- 
ing an ad hoc step in kinematic viscosity around r = 1 .8ro, using 
the same recipe as described in Masset et al. (2006), adopting a 
value of F - l.ol'^i - 10 for the ratio of the surface densi- 
ties outside and inside the cavity. We obtain a rather sharp IC of 
width Ar ^ O.Stq. Boundary conditions were chosen to be non- 
reflecting for the inner edge and such that there was a continuous 
outer source mass to maintain the surface density in the exter- 
nal regions of the disk (Masset et al. 2006 used non-reflecting 
boundary conditions for both edges). No planet was considered 
at this point. 

An initial one-dimensional run was performed to allow the 
cavity to form and reach a steady state configuration. The re- 
sulting density profile and total torque were analogous to those 
shown in Masset et al. (2006), showing the existence of a stable 
fixed point in the torque at a distance r* slightly larger than the 
nominal edge of the IC (Figure |2]l. A planet located outside r* 
would feel a negative torque and suffer a negative orbital decay 
towards the star Conversely, a body placed at a radial distance 
slightly smaller than r* would feel a positive torque and an out- 
ward radial migration. Thus, the point r = r* is a stable nesting 



place for the planet, as long as its mass is sufficiently low to 
avoid disruption of the density (and torque) profile. 

2.2. Test simulations of planet trapping 

After the disk density profile had stabilized, we included a planet 
of mass m in an initially circular orbit at r = 2.5ro. The system 
was allowed to evolve until the planet reached a stationary so- 
lution and the orbital decay effectively stopped. Figure [3] shows 
the results of two runs, the first with m - O.lMjup and the second 
with m = lOMjup. The location of the center of the IC is marked 
by the top horizontal dashed line, while the lower corresponds to 
the interior 2/1 mean-motion resonance (MMR) with the center 
of the IC. 

The smaller planet sufifers an almost constant orbital decay 
until it is halted at a location slightly outside the cavity, in accor- 
dance with the findings of Masset et al (2006). The more mas- 
sive planet, however, suffers a Type II migration that opens a 
gap in its co-orbital region and completely disrupts the density 
profile of the disk in its vicinity. The cavity edge is therefore 
incapable of generating a strong corotation torque and fails to 
trap the planet, which continues its orbital decay inside the cav- 
ity. The migration is finally stopped very close to the location of 
the 2/1 MMR with the IC where the diff'erential Lindblad torque 
reduces to zero, as predicted by Kuchner and Lecar (2002). 

In all cases, the eccentricities were damped, and the final or- 
bits were observed to be quasi-circular with e ~ 10"^. This final 
eccentricity is different from the results presented by Rice et al. 
(2008), where a significant eccentricity excitation is observed 
for Type II migration inside the inner disk edge. From what we 
have been able to see, the discrepancy arises from the different 
numerical setup. Rice et al (2008) truncate the disk at the mesh's 
inner edge, and the planets that they consider essentially orbit 
in a vacuum. This is in contrast to our situation in which some 
residual material surrounds the orbit. This material, by the ac- 
tion of co-orbital Lindblad resonances, is a powerful source of 
eccentricity damping. 

2.3. Scaling and comparison with observations 

We repeated the simulations for a total of 18 runs, with planetary 
masses ranging from 0.01 to 20 Jupiter masses, all starting from 
initial circular orbits at r = 2.5ro. The orange curve in the top 
plot of Figure |4] presents the final orbital period as function of 
the planetary mass. For comparison, the observational data are 
again drawn in circles. 

Since the location of our IC was in arbitrary units, we have 
a degree of freedom in the sense that we may shift the orange 
curve arbitrarily in the y-axis to fit the observed exoplanets. We 
find that the best fit gives an inner cavity edge located at approx- 
imately a ~ 0.03 AU corresponding to an orbital period of ~ 2 
days. This value is comparable with the inferred inner gas radii 
of T Tauri star disks from CO spectroscopy (Najita et al. 2007, 
Carr 2007) which seems to be located around 0.03 - 0.04 AU. 
Although our best-fit distance is slightly lower, any difference 
could be due to stellar parameters and the maximum CO veloc- 
ity (Carr 2007). 

Exoplanets with masses below one Jovian mass seem to have 
a lower-limit orbital period of about 2-3 days, consistent with 
the position of the equilibrium point being slightly beyond the 
IC (Masset et al. 2006). On the other hand, the orbital periods of 
more massive planets have a lower limit, this time in accordance 
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Fig. 4. Gray circles reproduce the mass-period distribution of 
known close-in exoplanets. The broad orange band represents 
the final semimajor obtained from the runs, scaled along the y- 
axis to fit the observational data. In the lower frame, a denotes 
the normalized semimajor axis defined in Eq. (|5]l. 



with the location of the 2/1 mean-motion resonance with the IC, 
as predicted by Kuchner & Lecar (2002). 



3. Tidal evolution of close-in exoplanets 

Although the results presented in the top graph of Figure |4] ap- 
pear to be encouraging, there are two important approximations 
that must be examined. First, the exoplanets orbit stars with dif- 
ferent radii and masses, and we have assumed that the scaling 
in the y-axis is the same. Although we have not assumed any 
origin for the inner cavity (MRI, stellar wind, etc.), it is almost 
certain that the location of the IC should be a function of the 
stellar mass. If we assume a very simple model in which r/c 
scales with then the position of both GJ1214b and GJ876d 
would be displaced above the broad orange curve, thus eliminat- 
ing their incongruity with respect to the simulations. However, 
CoRoT-7 has a stellar mass of 0.93Mq and its location in the 
diagram would still be conflicting. 

A second approximation is that we have neglected the 
later evolution of the planetary periods caused by stellar tides. 
Assuming almost circular planar orbits for the close-in planets, 
we can neglect the planetary tides and approximate the differen- 
tial equation for orbital decay as 



da 
dt 
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(see Ferraz-Mello et al. 2008, Jackson et al. 2009). This expres- 
sion can be easily integrated to yield 
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where gq - a(t - 0) is the initial value of the semimajor axis. 
After some simple algebraic manipulations, we can express this 
solution as 
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Equation O) is then independent of the stellar parameters, whose 
values are incorporated into the "normalized" semimajor axis a. 
We note that the rate of orbital decay is given by /3, which is 
linearly proportional to the planetary mass and inversely pro- 
portional to the stellar dissipation parameter Q'^. Is it of course 
probable that the tidal parameter itself depends on the stellar pa- 
rameters (see Barker and Ogilvie 2009), but there has so far been 
no clear indication of how it may vary nor by what magnitude. 

The lower plot in Figure |4] shows the distribution of a as 
a function of m for the close-in exoplanet population. The or- 
ange curve marks the position of the simulated lower limit to 
the semimajor axis. As expected, both GJ876d and GJ1214b are 
now above the lower limit, although CoRoT-7b still remains de- 
tached. We note that the bump in the distribution is still clearly 
visible. 

We can now use equation (O to check whether the present- 
day distribution could originate solely from tidal evolution. To 
test this idea, we generated a synthetic population of 10"* ficti- 
tious exoplanets distributed randomly in planetary mass across 
the interval m e [0.01, 12]Mjup and with random stellar ages T 
between one and eight Gyr. The normalized semimajor axis a of 
each planet was then evolved tidally for the age of the system 
using a stellar tidal parameter Q'^ - 10^. 

Figure|5]shows results considering three different initial dis- 
tributions, chosen to be uniform in orbital period P and log (m). 
In all cases, the lowest mass was equal to m = O.OlMjup, and 
the largest to ten Jovian masses. Since this preliminary analysis 
is only intended to be illustrative, we assumed solar-type stars 
for all bodies; consequently, the normalized semimajor axis a is 
equal to the nominal semimajor axis a. 

For the top plots, we assumed fictitious planets with orbital 
periods that have no lower limit. The initial data set is shown in 
the left plot, while the right-hand plot presents the final distribu- 
tion after evolution through stellar tides for T - S Gyr. Since the 
decay rate /3 is proportional to the planetary mass, more massive 
planets fall more rapidly towards the star, leading to a final pop- 
ulation that has a lack of large bodies with small orbital periods. 
This distribution shows little relation to the real distribution of 
planets. 

For the middle plots, we considered a lower limit to the or- 
bital periods equal to 3 days. This value is sufficiently large for 
terrestrial-type bodies to be virtually unaffected, although gi- 
ant planets still suffer significant orbital decay. Depending on 
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Fig. 5. Tidal evolution of fictitious planets (stellar tides alone) 
with - 10''. Initial distributions are shown on the left, while 
the evolved distributions are shown on the right. Stellar ages 
where chosen uniformly between one and eight Gyr 



their initial semimajor axis and stellar age, many massive plan- 
ets reach the Roche radius and are engulfed by the star, but a 
portion of the population remains. We note the final step-like dis- 
tribution, which is reminiscent of the real population. However, 
the bump is now located at ~ 0.3Mjup, thus at a lower mass than 
for the real exoplanets. 

Finally, in the lower plots of Figure |5] we considered an ini- 
tial population with a step in the orbital period, similar to that 
resulting from the hydrodynamical simulations. Even though the 
initial conditions are different from those depicted in the middle 
plots, there is no significant difference in the final distribution. 
This seems to indicate that any evidence of an initial disk-related 
structure would be smeared by the later tidal evolution. Thus, a 
bump in the present real population is not necessarily indicative 
of a dynamical structure prior to the dissipation of the gaseous 
disk. 

These results indicate the clear possibility that the observed 
bump in the exoplanet a-m distribution could be mainly due to 
stellar tidal effects. To test this proposal in more detail, we cal- 
culated how a lower bound in the values of a would be modified, 
for different values of Q^, after a timescale of T = 8 Gyr Results 
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Fig. 6. Distribution of close-in exoplanets in the a-m diagram. 
Gray curves show the end evolution due to stellar tides (at T = 8 
Gyr) of initial constant values of a. Each curve corresponds to 
a different ao- Stellar tidal parameters are given in the top right- 
hand side of each graph. The broad dashed curve in the top panel, 
corresponding to ao = 0.035, gives the best fit to the lower limit 
of the current exoplanet population. 



are shown in Figure|6]for three values of the stellar tidal param- 
eter Each gray curve shows the function 



.13/2 



2/13 



with r = 8Gyr 



(6) 



for different values of So- For comparison, each plot also repro- 
duces the present distribution of close-in exoplanets. 

The present-day lower limit to the normalized semimajor 
axis shows a very good agreement assuming ^ 0.035 and 
a tidal evolution with Q'* = 10^, represented by a bold dashed 
curve in the top graph. This value of the normalized semimajor 
axis corresponds to original orbital periods of P ~ 2.5 days for 
solar-type stars. We note that the other plots, corresponding to 
smaller values of do not show good correspondence for any 
adopted value of aq- This result appears to indicate that smaller 
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Fig. 7. Left: Current-day distribution of exoplanets in the a-m 
diagram. Black circles show planets for which values of M», 
and stellar ages T, are available. Gray circles show those bodies 
for which stellar data is missing. Right: Reconstructed origi- 
nal distribution assuming backward orbital evolution from stel- 
lar tides for the age of each parent star In both plots, the dashed 
curves show the lower limit to normalized semimajor axis, as a 
function of the planetary mass, obtained from the hydrodynami- 
cal simulations and scaled along the y-axis to fit the exoplanets. 
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Fig. 8. Left: Current values of the normalized semimajor axis 
a for close-in exoplanets as a function of the stellar age T,. 
Notice that bodies with a < 0.03 are limited to young systems 
(Tt < 3Gyr), while larger semimajor axes are found for all stellar 
ages. Riglit: Original distribution in the m-a diagram. Diagonal 
curves give limits to the orbital decay leading to disruption for 
two stellar ages: the continuous line show results for T, = 1 Gyr, 
while the dashed line corresponds to = 8 Gyr. 



values of Ql, and consequently faster orbital decays due to stel- 
lar tides, are inconsistent with the current population of close-in 
exoplanets. 

Although the timescale (T = 8Gyr) chosen for Figure |6]may 
seem arbitrary, equations (HJi and (|5]l show that the most relevant 
parameter for the orbital evolution is actually the ratio T/Q'^. 
Since the uncertainty in the tidal parameter is much larger than 
in the age of the system, it seems justified to use a single fixed 
value of T and to assume that no qualitative differences would 
be observed for other stellar ages. 

As a final test, we can attempt to reconstruct the original 
mass/semimajor-axis distribution of the exoplanet population, 
integrating each current value of a backwards in time using equa- 
tion (|6|l, adopting in this case the value of T equal to the age of 
each parent star T,. To perform this calculation, we require in 
addition to the stellar mass M, and radius estimates of each 
stellar age T,. However, this information is not available for all 
planetary systems. Out of the 133 original planets from our data 
set in direct orbits (i.e. eliminating the bodies believed to be in 
retrograde motion), we were only able to obtain a full set of stel- 
lar parameters for 94 exoplanets. Stellar data were obtained from 



the Simbad database for stellar properties and from Jackson et al. 
(2009). 

The present-day a-m distribution for this reduced population 
is shown as black circles in the left-hand plot of Figure |7] For 
comparison, gray circles show those planets for which complete 
stellar data is currently unavailable. The location of the lower 
bound to the semimajor axis deduced from the hydrosimulations 
is shown as a bold dashed curve. Unfortunately, many systems 
with giant planets with a < 0.025 do not have complete stellar 
properties, and the bump in the distribution is not clearly visible 
for the smaller population. 

The right-hand side plot of FigureQshows the "original" lo- 
cation of the smaller population after the backwards integration 
for each stellar age. For comparison, the inner-cavity-induced 
lower limit is again shown as a bold dashed curve, although 
the scaling in the y-axis has been modified to fit the new val- 
ues of semimajor axis. Although the distribution still shows a 
good agreement for planetary masses up to m ~ IMjup, most of 
the higher masses have increased their semimajor axis beyond 
a ~ 0.03, and little evidence remains of the bump. 

A possible explanation could be that most giant planets with 
small values of a were subsequently lost due to tidal disruption. 
To check this hypothesis, the left plot of Figure [8] shows the re- 
lation between the current values of the normalized semimajor 
axis as a function of the stellar age. Although planets with rel- 
atively large semimajor axes exist for all values of r», it seems 
very clear that giant planets with a < 0.025, which are primar- 
ily responsible for the bump in the distribution, belong to young 
systems with a maximum stellar age of J* ^ 3 Gyr. 

From Eq. (O, it is possible to estimate, for each planetary 
mass, the critical value of the normalized semimajor axis (i.e. 
Aq) that falls towards the star of a given stellar age T,. This is 
given approximately by 



2/13 



(7) 



which is simply obtained by setting the final semimajor axis to 
be equal to zero. Although the planet is believed to disrupt upon 
reaching the Roche radius and, thus, before impacting the star 
itself, as shown by Jackson et al. (2009) the orbital decay at such 
small semimajor axes is so swift that the timescales for both sce- 
narios are practically equal. 

The right-hand plot of Figure [8] once again reproduces the 
"original" distribution of exoplanets in the m-a diagram shown 
before in Figure|7] The diagonal lines show the values of Og, as a 
function of the planetary mass, for two stellar ages: - I Gyr 
(continuous line) and 7,-2, Gyr (dashed line). It is clear that 
even if giant planets were deposited by a hypothetical disk inner 
cavity at small values of the semimajor axis (a < 0.025), they 
would be rapidly absorbed by the parent star because of tidal 
effects, even for planetary systems as young as 1 Gyr. Thus, it 
is not unexpected that if any original bump in the distribution 
were created by disk-planet interactions, subsequent tidal effects 
would have eliminated most traces. 

4. CoRoT-7b 

CoRoT-7 appears to be a special case. As seen from the lower 
plot of Figure m its low mass and short orbital period mean that 
it is well separated from the m-a distribution observed for other 
close-in exoplanets. We recall that this planet orbits a solar-type 
star of mass M» ^ O.93M0. 

Jackson et al. (2009) proposed that these planetary bodies 
should be undergoing significant tidal evolution and orbital de- 
cay. Most bodies in this mass range would then be absorbed by 
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the star on timescales shorter than the age of the star According 
to this idea, CoRoT-7b owes its present existence solely to its 
star being relatively young (~ 1 .5 Gyr). 

However, numerical simulations of the tidal evolution of 
the CoRoT-7 planetary system (Ferraz-Mello et al. 2010) in- 
dicate that the current eccentricities should be extremely low 
(e «: 10"^), and that any primordial departure from circular 
motion would have been rapidly damped before any significant 
orbital decay occurred. Thus, tidal evolution would have been 
given primarily by stellar tides alone. In this case, the original 
location of CoRoT-7b in the primordial m-a distribution should 
be given by the right-hand plots of Figures [T] and [8j once again 
this planet appears to be detached from the remaining close-in 
planetary population. 

Since CoRoT7 harbors at least one additional planet, it is 
possible that mutual dynamical interactions might also explain 
this planet's proximity to the star. In multiple-planet systems, 
scattering is believed to have played an important role in sculpt- 
ing the general exoplanet distribution, and the same phenom- 
ena may have occurred in this system. According to this idea, 
CoRoT7-b could have had a close encounter with CoRoT7- 
c (or with an additional ejected planet) and suffered a signif- 
icant reduction in its semimajor axis. Subsequent tidal inter- 
actions would have circularized its orbit to its present state. 
However, CoRoT7 is not the only multiple-planet system in this 
region. Both the GI581 and HD40307 planetary systems have 
two known super-Earths with short orbital periods that never- 
theless lie above the expected lower limit. Hence once again 
CoRoT7 seems to be different. 

A possible explanation may lie elsewhere. Valencia et al. 
(2010) and Jackson et al. (2010) proposed that CoRoT-7b could 
be the solid core of a primordial giant planet whose gaseous en- 
velope was lost due to evaporation. A recent re-evaluation of the 
radial velocity data by several authors indicate that the mass of 
CoRoT-7b could be as high as 9Me (Ferraz-Mello et al. 2010) 
or as low as 2M^ (Pont et al. 2010). Although a high current 
mass appears indicative of a rocky/iron composition and that the 
original mass was at most twice the present value (Valencia et 
al. 2010), a lower mass is consistent with a lighter composition 
and a envelope-depleted gas giant. 

If CoRoT-7b were indeed the solid core of a primordial gas 
giant, then its location in the m-a diagram (right-hand plot of 
Figure |7]i would be located close to the dashed curve and, thus, 
consistent with the rest of the planetary population. We still need 
of course to explain how this planet suffered an evaporation of its 
gaseous envelope on timescales shorter than one Gyr, and why it 
appears to be the only example of this effect. 



5. Conclusions 

We have attempted to understand the dynamical origin and evo- 
lution of the mass-period distribution of close-in exoplanets. The 
present-day population shows a distinctive discontinuity located 
at approximately one Jovian mass. Smaller planets have orbital 
periods longer than P ~ 2.5 days, while higher masses are found 
to have periods as short as f ~ 1 day. 

We have found that the combined effects of tidal evolution 
and disk-planet interactions with an inner cavity (IC) in the 
gaseous disk can explain most of the observed characteristics. 
The current distribution appears to be compatible with an in- 
ner disk edge located approximately at distances of a ^ 0.035, 
which for solar-type stars corresponds roughly to orbital periods 
of P ^ 2.5 days. This value is consistent with the inner gas radii 



for T Tauri stars as estimated from CO spectroscopy (Najita et 
al. 2007, Can- 2007). 

Planets below a certain critical mass ~ Mj^p are trapped 
just outside the IC as found by Masset et al. (2006). The loca- 
tion of the stationary solution with respect to the IC is practi- 
cally mass-independent. In contrast, bodies with m > nic enter 
a regime characterized by a Type II migration that causes sig- 
nificant perturbations to the density profile of the disk; conse- 
quently the IC edge cannot generate a significant positive coro- 
tational torque and does not stop the orbital decay. As predicted 
by Kuchner and Lecar (2002), migration only brakes inside the 
inner edge at a 2/ 1 mean-motion resonance with the cavity edge. 

For reasonable values of the disk viscosity, we expect a gap 
opening to occur when the height of the disk is approximately 
equal to the Hill radius of the planet. For a scale height equal 
to HI r ^ 0.05, this implies a minimum mass of m ^ 0.4Mjup, 
a value similar to our critical mass m^. Adopting a value of 
one Jovian mass for the critical mass leads to a slightly larger 
value of HI r ^ 0.07. However, given the uncertainties involved 
in both the gap opening criteria and its dependence on the vis- 
cosity, the values may be considered to be virtually equivalent. 
Thus, the location of the bump in the observed distribution of 
close-in planets appears to be consistent with a mass thresh- 
old for gap opening (Crida et al. 2006 and references therein). 
Although the present data are sparse and plagued by the ad- 
ditional effects described in this work, we may expect that the 
aforementioned results will be confirmed by further detections. 
In contrast, when the statistics become sufficiently robust, the 
location of this bump may be used to place constraints on the 
physical properties of the inner disk, in terms of both tempera- 
ture and effective viscosity, as these quantities feature in the gap 
opening criterion. 

Since the edge of the IC created in our simulations is placed 
at an arbitrary distance from the star, we have a certain degree of 
freedom when fitting the numerical mass-period distribution to 
the real planets (as seen in Figure|4]l. We have adopted a vertical 
displacement that minimizes the number of exoplanet inside the 
cavity edge, but this is not the only option. It may be argued that 
it would be better to fit the synthetic curve with the location of 
sub-giants (i.e. m ~ 0.5Mjup) in the m-P diagram. We note that 
these bodies exhibit a smaller dispersion in orbital period than 
observed for any other mass range, and can be seen as a com- 
pact group in both plots of Figure|4] However, this displacement 
would lead to our accepting a larger number of small planets 
inside the cavity edge, bodies whose subsequent tidal evolution 
should have been very small. 

Whatever the choice, the qualitative results are not affected. 
Moreover, the ratio of the stopping values of P (for small and 
large bodies) is scale-independent, and found to be slightly 
larger than 2. This is because although higher masses are stopped 
in a 2/1 MMR with the IC, the small bodies are trapped outside 
the disk edge. Once again, the distribution of real planets seems 
to yield a similar ratio. 

The subsequent tidal evolution of the close-in planets is con- 
sistent with a stellar tidal parameter of Q'^ - 10^, a value simi- 
lar to that predicted by Schlaufman et al. (2010) from synthetic 
population models. Smaller parameters, leading to higher rates 
of orbital decay, do not lead to distributions similar to the ob- 
served population. This is also consistent with the analysis of 
Ogilvie & Lin (2007). A consequence of the tidal evolution is 
the removal of most of the original gas giants with short orbital 
periods and their substitution by exoplanets that were originally 
farther away. Thus, we expect that many of the primordial plan- 
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ets with P < 2 days and m > IMjup might have been tidally 
disrapted and absorbed by their parent stars. 

Although this scenario is consistent with the properties of 
most of the exoplanetary population, it appears difficult to ex- 
plain the present-day mass and orbit of CoRoT-7b. A possible 
explanation is to assume that the planet started its life as a gas 
giant whose gas envelope was completely evaporated (Valencia 
et al. 2010, Jackson et al. 2010). 

Last of all, in the present scenario we have neglected the role 
of the orbital inclinations; however, the same results should be 
expected as long as the inclinations are not very large. Three- 
dimensional studies of disk-planet interactions in the linear ap- 
proximation (e.g. Tanaka and Ward 2004), as well as numerical 
simulations (e.g. Cresswell et al. 2007), show little effect on a 
finite inclination on the migration timescale. Similar results are 
also found for orbital evolution due to tidal eff'ects (e.g. Ferraz- 
Mello et al. 2008, Barker and Ogilvie 2009). For planets in ret- 
rograde orbits results are, however, is more difficult to evaluate. 
It is unclear whether gas disks in retrograde motion (with re- 
spect to the stellar spin) would have inner cavities or how such 
a structure would interact with planets in its vicinity. Similarly, 
the tidal evolution of bodies in retrograde motion is poorly un- 
derstood, thus it is impossible at present to ascertain how the 
results of this work could be extended to these systems. 
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